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Abstract 


Defence  Research  and  Development  Canada  iValcarticr  (DRDC  Valcartier)  has  been  developing 
protection  systems  for  many  years  to  reduce  the  vulnerability  of  light  armoured  vehicles  (LAV) 
against  mine  blast.  To  assist  in  the  development  of  these  protections,  experimental  tests  and  finite 
element  (FE)  analyses  are  performed.  To  carry  out  these  numerical  calculations,  initial  conditions 
such  as  the  loading  prescribed  by  a  mine  on  a  structure  need  to  be  simulated  adequately.  The 
effects  of  blast  on  structures  depend  often  on  how  these  initial  conditions  are  estimated  and 
applied.  In  this  report,  two  methods  were  used  to  simulate  a  mine  blast:  the  arbitrary  Lagrangian- 
Eulerian  (ALE)  and  the  smoothed  particle  hydrodynamics  (SPH)  formulations.  The  comparative 
study  was  done  for  a  simple  and  a  more  complex  target. 

The  first  target  is  an  aluminum  plate  placed  on  four  steel  legs  and  centred  over  a  surrogate  mine 
filled  with  6  kg  of  C4  which  represents  a  typical  buried  mine  blast  scenario.  Two  FE  models 
were  generated:  an  ALE  model  and  a  SPH  model.  The  final  deformation  of  the  plate  was 
measured  for  both  approaches  and  was  compared  with  experimental  measurements. 

The  second  target  is  a  mock-up,  representing  a  section  of  the  side  of  a  typical  LAV,  subjected  to  a 
typical  improvised  explosive  device  (IED).  Two  FE  models  were  generated:  an  ALE  model  and  a 
SPH  model.  Parametric  studies  were  done  on  both  models  and  the  best  results  were  compared  to 
the  experimental  ones.  Each  comparison  included  the  velocity  at  the  center  of  the  sponson  wall, 
sponson  top  and  sponson  sidewall. 
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R  sum 


Le  centre  Recherche  et  d  veloppemcnt  pour  la  defense  Canada  PI  Val carticr  (RDDC  Valcartier) 
dfkeloppe  des  systQnes  de  protection  depuis  plusieurs  annas  pour  rQluire  la  vulnQabilitndes 
vfhicules  blinds  lexers  (VBL)  contre  les  mines  Deffet  de  souffle.  Pour  aider  au  dfkeloppement 
de  ces  protections,  des  tests  exp  rimentaux  et  des  analyses  par  1  merits  finis  (EF)  sont  r  alisS. 
Pour  effectuer  les  simulations  num  Piques,  des  conditions  initiales  telles  que  lSffet  de  souffle  de 
lSxplosion  d  une  mine  sur  une  structure  doivent  tre  simulfes  ad  quatement.  LSffet  de  souffle 
de  lSxplosion  dfjrend  souvent  de  la  faQn  dont  les  conditions  initiales  sont  esti ruffes  et 
appliqu  es.  Dans  ce  rapport,  deux  m  thodcs  ont  1  util  is  es  pour  simuler  lSffet  de  souffle  de 
lSxplosion  d  une  mine  :  la  formulation  Lagrangienne  Eul  rienne  arbitraire  (ALE)  et  la 
formulation  hydrodynamique  des  particules  lisses  (SPH).  LStude  comparative  a  IDrl  at  is  e  pour 
une  cible  simple  et  une  cible  plus  complexe. 

La  premise  cible  est  une  plaque  d aluminium  placffe  sur  quatre  poutres  dScier  et  centra  au- 
dessus  d  une  mine  Simula  de  6  kg  de  C4,  ce  qui  repr Sente  un  sc  nario  typique  de  mine  enfouie. 
Deux  modfles  EF  ont  3f]gfhQS  :  un  modfle  ALE  et  un  modfle  SPH.  La  deformation  finale  de 
la  plaque  a  3D  mesura  pour  les  deux  approches  et  elle  a  □□  compara  aux  mesures 
exp  Qimentales. 

La  deuxi  me  cible  est  un  mod  le  partiel  repr  Scntant  une  section  d  un  c  tHdun  VBL  typique, 
soumis  □  lSffet  de  souffle  d  un  engin  explosif  improvise (EEI)  typique.  Deux  mod  les  EF  ont 
□□  gfhQS  :  un  modQe  ALE  et  un  modQe  SPH.  Des  3udes  paramftriques  ont  3GJ 
effectuSs  pour  les  deux  mod  les  et  les  meilleurs  r  sultats  ont  □□  comparS  aux  r  sultats 
exp  Qimentaux.  Pour  chacune  de  ses  3udes,  la  vitesse  au  milieu  du  dessus  de  Mile,  du  mur  au- 
dessus  de  Mile  et  du  mur  en-dessous  de  Mile  du  modfle  partiel  ont  aDcomparas  aux  mesures 
exp  Qimentales. 
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Executive  summary 


Comparison  of  ALE  and  SPH  methods  for  simulating  mine  blast 
effects  on  structures 

Toussaint  G.,  Bouamoul  A.;  DRDC  Valcartier  TR  2010-326;  Defence  R&D 

Canada  QValcartier;  December  2010. 

Background:  Light  armoured  vehicles  (LAV)  are  exposed  to  non  conventional  threats  such  as 
improvised  explosive  devices  (IED).  To  improve  the  crew  survivability,  armoured  vehicles 
structures  need  to  be  reassessed.  To  achieve  this  goal,  experimental  setups  are  developed  and 
experimental  testing  is  conducted  at  the  Defence  Research  and  Development  Canada  □Valcartier 
(DRDC  Valcartier)  center  to  evaluate  armoured  vehicle  structure  subjected  to  blast  loadings.  In 
conjunction  with  these  tests,  finite  element  (FE)  models  are  used  to  optimise  the  design  of 
protection  systems.  However,  to  perform  numerical  calculations,  initial  conditions  such  as  the 
loading  prescribed  by  a  mine  on  a  structure  need  to  be  simulated  adequately.  The  effects  of  blast 
on  structures  depend  often  on  how  these  initial  conditions  are  estimated  and  applied. 

Different  approaches  such  as:  CONWEP,  Westine,  CHINOOK  and  mine  pressure  model  were 
investigated  by  DRDC  Valcartier  to  simulate  mine  blast  on  structures.  The  intent  of  this  report  is 
to  continue  this  investigation  by  evaluating  two  approaches  available  in  LS-DYNA  code:  the 
arbitrary  Lagrangian-Eulerian  (ALE)  and  the  smoothed  particle  hydrodynamics  (SPH)  methods. 
Two  scenarios  were  simulated:  a  target  plate  and  a  mock-up  structure  (representing  the  side  of  a 
typical  LAV).  The  numerical  results  obtained  in  both  scenarios  were  compared  to  the 
experimental  results. 

Principal  results:  For  the  plate  setup,  both  approaches,  ALE  and  SPH,  were  suitable  to  simulate 
the  effect  of  a  buried  mine.  However,  SPH  is  the  approach  that  was  in  best  accordance  with  the 
final  experimental  deformation.  In  the  mock-up  scenario,  the  ALE  model  better  reflected  the 
evolution  of  the  air  blast  than  the  SPH  model  did.  The  relevant  aspects  to  consider  when  using 
the  ALE  method  to  model  air  blast  are,  to  use  a  fixed  domain  with  a  good  mesh  and  to  choose 
judiciously  the  air  material  properties  and  its  corresponding  equation  of  state.  In  general,  the 
factors  that  may  influence  an  ALE  or  SPH  results  are  many.  Among  them  are  the  complexity  and 
dimensions  of  the  structures  under  the  blast,  the  mesh,  the  material  constitutive  models  and  their 
corresponding  equation  of  state  and  the  boundary  conditions.  All  these  parameters  need  to  be 
studied  in  details  when  modeling  blast. 

Significance  of  results:  The  results  obtained  in  this  report  will  help  to  better  simulate  the  mine 
blast  effects  on  different  structures  using  the  ALE  and  SPH  methods  and  consequently,  it  will 
improve  the  analysis  of  protection  systems  using  FE  method. 


iii 


DRDC  Valcartier  TR  2010-326 


Sommaire 


Comparison  of  ALE  and  SPH  methods  for  simulating  mine  blast 
effects  on  structures 

Toussaint  G.,  Bouamoul  A.;  DRDC  Valcartier  TR  2010-326;  R  et  D  pour  la 

defense  Canada  QValcartier;  DBcembre  2010. 

Contexte  :  Les  v  hiculcs  blind  s  I  gcrs  (VBL)  sont  exposV  Ddes  menaces  non  conventionnelles, 
telles  que  les  engins  explosifs  improvise  (EEI).  Dans  le  but  d  augmcntcr  la  survie  des 
occupants,  la  structure  des  v  hiculcs  blinds  doit  trc  r  Bxamin  t.  Pour  atteindre  ce  but,  des 
montages  exp Bimentaux  sont  d  veloppB  et  des  tests  exp  Bimentaux  sont  r  alisV  au  centre 
Recherche  et  d  vcioppcmcnt  pour  la  defense  Canada  □  Valcartier  (RDDC  Valcartier)  pour 
[Valuer  la  structure  des  v  hiculcs  blind  s  soumis  au  souffle  de  lBxplosion.  En  m  me  temps,  des 
mod  des  d  OH  Orients  finis  (EF)  sont  util  i  s  s  pour  optimiser  le  design  de  syst  mes  de  protection. 
Toutefois,  pour  r  aliscr  les  calculs  num  riques,  des  conditions  initiales  telles  que  lBffet  de  souffle 
de  lBxplosion  d  une  mine  sur  une  structure  doivent  Ore  simul  es  ad  quatement.  LBffet  de 
souffle  de  lBxplosion  defend  souvent  de  la  fa  Bin  dont  les  conditions  initiales  sont  estimBs  et 
appliquBs. 

Diff  rentes  approches  telles  que  :  CONWEP,  Westine,  CEIINOOK  et  le  mod  le  de  pression  ont 
!□  tudi  s  par  RDDC  Valcartier  pour  simuler  lBffet  de  souffle  de  lBxplosion  sur  les  structures. 
LBbjectif  de  ce  rapport  est  de  poursuivre  cette  Bude  en  Valuant  deux  mBhodes  disponibles  dans 
le  code  LS-DYNA  :  la  formulation  Lagrangienne  Eul Bienne  arbitraire  (ALE)  et  la  formulation 
hydrodynamique  des  particules  lisses  (SPEl)  pour  simuler  lBffet  de  souffle  de  lBxplosion  d  une 
mine  sur  une  structure  typique.  Deux  scBiarios  ont  BDsimulV  :  une  plaque  cible  et  un  modBe 
partiel  d  un  v  hiculc  blind  (repr  sentant  une  section  d  un  cBHd  un  VBL  typique).  Les  r  sultats 
numBiques  obtenus  dans  chaque  scBiario  ont  1  icompar  s  aux  r  sultats  exp Bimentaux. 

Principaux  rVultats  :  Dans  le  cas  de  la  plaque  cible,  les  deux  approches,  ALE  et  SPH,  peuvent 
Bre  utilisBs  pour  simuler  lBffet  de  souffle  d  une  mine  enfouie.  Toutefois,  SPH  est  Mpproche 
dont  le  rl  sultat  Bait  le  plus  conforme  avec  la  deformation  expBimentale  finale.  Pour  le  modBe 
partiel,  le  modBe  ALE  a  mieux  reflBD involution  du  souffle  de  lBxplosion  dans  Mir  que  le 
modBe  SPH.  Les  aspects  les  plus  importants  dont  il  faut  tenir  compte  lorsque  le  modBe  ALE  est 
utilisDpour  modBiser  un  souffle  de  lBxplosion  dans  Mir  sont  d  utiliscr  un  domaine  fixe  avec  un 
bon  maillage  et  de  choisir  judicieusement  les  propriBV  des  matBiaux  et  lCHpiation  dffiat 
correspondante.  En  gBiBal,  plusieurs  facteurs  peuvent  influencer  les  r  Sultats  ALE  et  SPH, 
notamment,  la  complcxitlet  les  dimensions  des  structures  soumises  au  souffle  de  lBxplosion,  le 
maillage,  les  modBes  constitutifs  des  matBiaux  et  leur  quation  d  'tat  correspondante  et  les 
conditions  aux  frontiBes.  Tous  ces  paramBres  doivent  Bre  tudi  &  en  dBail  lorsquBn  simule 
lBffet  de  souffle  de  lBxplosion. 

Signification  des  r  sultats  :  Les  r  sultats  obtenus  dans  ce  rapport  vont  aider  B mieux  simuler 
lBffet  de  souffle  de  lBxplosion  d  une  mine  sur  diffBentes  structures  en  utilisant  les  mBhodes 
ALE  et  SPH  et,  par  consBjuent,  cela  amBiorera  Mnalyse  des  systHnes  de  protection  par  la 
mBhode  des  EF. 
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1  Introduction 


Defence  Research  and  Development  Canada  iValcarticr  (DRDC  Valcartier)  has  been  developing 
protection  systems  for  many  years  to  reduce  the  vulnerability  of  light  armoured  vehicles  (LAV)  to 
mine  blast.  To  assist  in  the  development  of  these  protection  systems,  DRDC  Valcartier  has 
developed  experimental  tests  and  has  acquired  an  expertise  in  finite  element  (FE)  analyses 
\Toussaint  and  Durocher,  2008,  Bouamoul  and  Gaudreault,  2006,  Durocher  et  al.,  2006, 
Bouamoul  and  Durocher,  2005,  Williams  and  Fillion-Gourdeau,  2002,  Williams  et  al.,  2002, 
Williams  and  Boon,  2000,  Williams,  1999].  To  perform  finite  element  calculations,  initial 
conditions  such  as  the  prescribed  loading  generated  by  a  mine  on  a  vehicle  structure  need  to  be 
simulated  adequately.  Accurately  simulating  the  effects  of  blast  on  structures  often  depends  on 
how  these  initial  conditions  are  estimated  and  applied. 

To  simulate  the  loading  conditions  generated  by  mine  blast,  different  approaches  were  explored 
in  recent  years.  This  report  briefly  describes  these  approaches,  which  include  the  CONWEP 
model,  the  Westine  model,  the  pressure-based  mine  loading  model,  the  CHINOOK  code,  the 
arbitrary  Lagrangian-Eulerian  (ALE)  method  and  the  smoothed  particle  hydrodynamics  (SPH) 
method.  Because  each  of  these  techniques  has  its  own  limitations,  DRDC  Valcartier  continues  to 
investigate  and  develop  new  features  to  predict  mine  blast  and  to  model  the  interaction  of  blast 
with  different  structures. 

The  aim  of  this  report  is  to  compare  numerical  results  obtained  using  the  two  FE  formulations 
ALE  and  SPH  with  experimental  data.  The  comparison  was  done  for  two  scenarios  and  the 
results  obtained  numerically  were  compared  with  the  experimental  data.  The  first  scenario, 
representing  a  typical  buried  mine  blast,  had  an  aluminum  plate  placed  on  four  steel  legs  and 
centred  over  a  surrogate  mine  filled  with  6  kg  of  C4.  The  experimental  setup,  the  FE  models  and 
the  results  are  all  given  in  Chapter  2.  The  second  scenario,  a  mock-up  representing  a  section  of  a 
typical  LAV,  was  subjected  to  a  typical  improvised  explosive  device  (IED)  air  blast.  The 
experimental  setup,  the  FE  model  and  the  results  are  detailed  in  Chapter  3.  For  this  complex 
scenario,  several  parametric  studies  were  conducted  for  the  ALE  and  SPH  methods.  For  the  ALE 
model,  mesh  sensitivity  analysis  of  the  mock-up  and  air  were  conducted.  A  comparison  between 
a  moving  and  a  fixed  air  domain  was  done  and  the  effect  of  the  explosive  geometry  was  studied. 
The  effect  of  the  air  material  properties  and  equation  of  state  were  also  investigated.  For  the  SPH 
model,  a  mesh  sensitivity  analysis  was  conducted  and  the  results  obtained  with  several 
formulations  of  the  particle  approximation  theory  were  compared.  Finally,  using  the  best 
combination  of  parameters  obtained  for  each  models,  a  comparison  was  made  between  the  ALE 
and  the  SPH  finite  element  results  and  the  experimental  data. 

This  work  was  done  at  DRDC  Valcartier  between  June  2007  and  September  2010,  under  the 
project  12ro03. 
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1.1  CONWEP  model 


CONWEP  is  a  collection  of  conventional  weapons  effects  calculations  from  the  equations  and 
curves  of  TM  5-855-1  [Randers-Pehrson  and  Bannister,  1997],  In  LS-DYNA  [LS-DYNA 
Keyword  User H  Manual,  2007 ]  FE  code,  the  *Load_blast  feature  was  implemented  based  on 
these  equations  and  curves  [Kingery  and  Bulmash,  1984 ]  and  was  used  at  DRDC  Valcartier  by 
Williams  et  al.  [2002].  CONWEP  model  is  valid  for  free  air  detonation  of  a  spherical  charge  and 
for  a  surface  detonation  of  a  hemispherical  charge.  However,  it  has  some  limitations  since  it 
considers  only  the  explosive  mass  but  not  soil  moisture  and  type. 

According  to  Williams  et  al.  [2002],  when  compared  with  experiments,  the  CONWEP  model 
showed  poor  results  that  could  be  explained  by  the  confining  effect  of  the  soil.  However,  the 
weight  of  the  charge  was  scaled  to  fit  the  experiment  and  should  be  rescaled  for  each  new  charge. 
According  to  Slavik  [2009-A],  shadowing  and  focusing  of  the  blast  are  not  taken  into  account  by 
this  model. 


1.2  Westine  model 

The  development  of  the  Westine  model  was  based  on  two  charges:  0.0567  kg  (1/8  lb)  and  0.2268 
kg  (1/2  lb)  [  Westine  et  al.,  1985].  To  use  this  empirical  equation,  four  criteria  which  depend  on 
the  mine  depth  of  burial,  the  standoff  distance  from  the  centre  of  the  mine,  the  energy  release  in 
the  explosive  charge,  the  cross-sectional  area  of  the  mine,  the  lateral  distance  to  the  centre  of  the 
mine,  the  soil  density  and  the  seismic  P-wave  velocity  in  the  soil,  must  be  satisfied  [  Westine  et 
al.,  1985,  Tremblay,  1998].  For  typical  land  mine  charges  (e.g.  6  kg  to  10  kg),  the  Westine  model 
was  extrapolated  linearly  with  respect  to  the  weight  of  the  mine  [  Williams  et  al.,  2002]. 

The  US  Army  TACOM  Impulse  model  developed  by  Westine  et  al.  [1985]  is  an  improvement 
from  the  CONWEP  model,  since  some  parameters  like  the  mine  depth  of  burial,  charge  size, 
target  height,  and  soil  density  are  taken  into  account. 

The  Westine  empirical  equation  was  incorporated  in  MinePre,  a  software  under  development  at 
DRDC  Valcartier  for  mine  blast  loading.  The  diameter  of  mine,  mine  mass,  TNT  equivalence, 
location  of  the  centre  of  the  mine,  depth  of  burial  and  soil  density  are  the  parameters  that  can  be 
used  as  an  input  to  compute  the  loading  profile.  The  soil  conditions  and  type  are  not  taken  into 
account  by  this  model. 

Williams  et  al.  [2002]  performed  numerical  simulations  on  blast  mine  effects  on  armoured 
vehicles  using  the  Westine  model  and  showed  that  the  impulse  loading  computed  by  the  model 
needed  to  be  scaled  to  fit  the  experimental  results.  The  authors  concluded  that  it  was  probably 
due  to  some  soil  conditions  that  are  not  taken  into  account  by  the  model. 

The  Westine  model  is  useful  and  accurate  for  loading  a  simple  geometry  like  flat  plate.  However, 
in  cases  where  the  geometry  is  complex  or  when  complex  interactions  between  the  structure  and 
the  soil  must  be  taken  into  account,  approaches  other  than  the  Westine  model  need  to  be 
considered. 


2 


DRDC  Valcartier  TR  2010-326 


1.3  Pressure-based  mine  loading  model 


The  empirical  equation  developed  by  Westine  et  al.  [1985]  was  used  to  predict  the  impulsive 
loading  caused  by  a  landmine  detonation.  One  limitation  of  this  impulse  model  is  its  inability  to 
load  a  FE  structure  made  of  solid  elements.  The  pressure-based  mine  loading  model  was 
developed  to  address  this  limitation.  The  pressure  model,  which  is  mainly  based  on  the  Westine 
impulse  model,  uses  a  pressure  time-space  distribution  to  apply  the  loading  conditions  on  the 
structure.  A  series  of  numerical  simulations  were  conducted  to  define  the  appropriate  pressure¬ 
time  history  to  be  applied  to  elements.  The  model  was  also  incorporated  in  the  MinePre  software. 

Due  to  the  limitations  of  the  last  three  models,  other  methods  like  the  CFIINOOK  code,  the  ALE 
and  the  SPE1  formulations  were  investigated. 

1.4  CHINOOK  code 

CFIINOOK  is  a  computational  fluid  dynamics  (CFD)  code  developed  by  Martec  Limited 
www.martec.com1.  It  has  the  ability  to  model  soil  effects,  which  is  very  important  in  landmine 
blast  applications,  and  includes  multi-phase  material  models.  Previous  studies  were  performed  to 
compare  this  CFD  model  against  experimental  results  [ Donahue  et  al.,  2005,  Donahue  et  al., 
2004,  Donahue  et  al.,  2003]  and  also  for  design  puipose  [Martec  Limited,  2006].  In  general,  the 
CHINOOK  code  has  the  ability  to  capture  some  of  the  physics  involved  in  the  simulation  such  as 
the  effects  of  blast  clearing,  focusing,  sheltering  and  channelling,  which  simple  codes  usually 
cannot  do.  However,  the  structure  is  usually  modelled  as  a  rigid  body  and  this  code  takes  more 
time  to  run  than  simple  codes. 

1.5  ALE  method 

An  alternative  to  the  Lagrangian  formulations  is  to  use  an  arbitrary  Lagrangian  Eulerian  (ALE) 
formulation  to  simulate  a  mine  detonation.  The  ALE  formulation  has  the  ability  to  resolve 
problems  with  large  mesh  distortion  [Donea  et  al.,  2004].  In  this  formulation,  the  reference 
domain  is  moving  and,  a  particular  case  is  the  Eulerian  formulation  where  the  reference  domain  is 
fixed.  The  ALE  method  is  divided  into  two  steps.  In  the  first  phase,  the  mesh  is  deformed  as  if  it 
was  Lagrangian  and  in  the  second  phase,  called  the  advection  phase,  computations  at  the  element 
boundaries  are  performed.  It  is  like  remapping  the  deformed  mesh  to  the  original  fixed  mesh 
[Alia  and  Souli,  2006]. 

The  literature  provides  different  examples  of  simulation  of  an  explosion  using  the  ALE  method 
using  LS-DYNA  code  [CendOi  et  al.,  2007,  Li  et  al.,  2007,  Larsen  and  Jorgensen,  2007,  Alia  and 
Souli,  2006,  Mahmadi  et  al.,  2004,  Mullin  and  O  Toole,  2004,  Vulitsk y  and  Kami,  2002,  Hilding, 
2001,  Wang,  2001].  The  ALE  method  allows  for  shadowing  and  focusing  of  a  fluid  (e.g.  blast, 
air);  however  it  has  some  limitations.  For  example,  to  accurately  capture  the  shock  wave  profile, 
fine  mesh  must  be  used;  therefore  computing  time  may  be  very  long  when  large  structures  are 
involved  or  when  the  explosive  is  located  a  few  metres  from  the  structure.  Another  limitation  is 
that  there  are  many  parameters  in  the  LS-DYNA  code  to  work  with.  It  is  also  possible  to 
experience  advection  errors  in  the  model,  which  results  in  leakage  into  the  structure  [LS-DYNA 
Keyword  Users  Manual  Version  971,  2007]. 
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Recently,  the  *Load _blast  enhanced  feature  [ Slavik ,  2009-A,  Slavik,  2009-B  and  Morten,  2008 ] 
was  added  to  the  LS-DYNA  code.  One  option  consists  in  applying  the  blast  pressure  profile 
generated  by  the  feature  to  the  closed  ALE  air  that  surrounds  the  target.  It  is  advantageous  to  use 
this  capability  when  the  blast  source  is  far  from  the  target  because,  while  still  addressing 
shadowing  and  focusing  of  the  blast,  it  is  not  necessary  to  model  the  complete  air  domain  from 
the  blast  source  to  the  target,  and  therefore,  the  computation  time  is  reduced  significantly. 
Further  numerical  simulations  are  required  to  verify  the  limitations  of  this  method. 

1.6  SPH  method 

The  smoothed  particle  hydrodynamics  (SPH)  method  was  initially  developed  by  Lucy  [1977], 
Gingold  and  Monaghan  [1977]  to  solve  problems  in  astrophysics  [Bala,  2007].  Lacome  [2000] 
has  implemented  this  method  in  LS-DYNA  code  to  avoid  problems  associated  with  large  mesh 
deformations  and  tangling  that  occurred  in  high  velocity  impact  problems.  In  SPH  method,  the 
mesh  is  represented  by  a  set  of  particles  moving  like  a  fluid.  To  provide  the  solution  of  the  entire 
problem,  governing  equations  are  resolved  for  each  particle  [Lacome,  2008,  Lacome,  2000,  LS- 
DYNA  Theory  Manual,  2006]. 

The  SPH  approach  was  used  to  predict  the  structural  response  of  a  plate  to  a  land  mine  detonation 
[Lacome,  2008]  and  to  predict  the  structural  response  of  a  mock-up  of  a  typical  LAV  to  an  air 
blast  [Toussaint  and  Durocher,  2008].  In  these  two  studies,  close  correlation  with  the 
experiments  was  achieved  with  the  SPH  approach. 
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2  First  scenario:  a  target  plate 


2.1  Experimental  setup 

In  order  to  validate  the  ALE  and  SPH  approaches,  simulations  of  an  experiment  performed  at 
DRDC  Valcartier,  using  the  experimental  blast  plate  configuration,  shown  in  Figure  1,  were 
performed.  A  1.828  m  x  1.828  m  x  0.03175  m  (72  in  x  72  in  x  1.25  in)  aluminum  5083-H131 
plate  was  placed  on  four  steel  legs  and  centred  over  a  surrogate  mine  containing  6  kg  of  C4 
[Donahue  et  al.,  2005] .  Lengths  of  steel  tubing  having  a  1  ft  by  1  ft  section  was  placed  around 
the  perimeter  of  the  plate.  Mass  was  added  to  the  top  of  the  steel  tubing  to  anchor  the  apparatus. 
The  mine  was  buried  in  a  compacted-gravel  roadbed  with  a  moisture  content  of  5%  and  density  of 
2300  kg/m3.  The  landmine  was  placed  with  an  overburden  of  5  cm  and  the  pressure  applied  to 
the  plate  at  a  standoff  distance  of  41  cm  was  measured.  This  setup  was  ideal  for  a  comparison 
between  the  ALE  and  SPH  methods  because  the  flat  plate  loading  area  and  the  simple  target 
geometry  would  not  favour  one  approach  over  the  other. 


Figure  1:  DRDC  experimental  plate  set-up 


2.2  Finite  element  models 

The  LS-DYNA  code  was  used  to  model  the  structural  response  of  the  aluminum  plate  to  6  kg  of 
C4.  To  reduce  the  number  of  elements  and  nodes  and  hence  the  calculation  time,  only  a  quarter 
of  the  apparatus,  as  shown  in  Figure  2,  was  modelled.  This  was  done  in  order  to  take  advantage 
of  the  planes  of  symmetry.  The  ALE  and  SPH  models  of  the  test  apparatus  developed  for  this 
calculation  are  given  in  Figure  2. 

In  both  ALE  and  SPH  models,  the  target  (plate)  and  the  steel  frame  were  meshed  respectively 
with  1,296  and  4,024  shell  elements.  Belytschko-Tsay  shell  elements  were  used  with  five 
interpolation  points  through  the  thickness.  Lumped  element  masses  were  used  to  model  the 
added  mass  on  top  of  the  steel  tubing. 

In  the  ALE  model,  the  air,  soil  and  mine  were  modelled  respectively  with  82,524,  56,936  and  196 
hexahedron  ALE  elements.  In  the  SPH  model,  the  soil  and  the  mine  were  modelled  with  196,010 
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and  8,000  SPH  elements  respectively  and  the  default  particle  approximation  was  used.  A  refined 
mesh  was  defined  for  the  region  of  the  soil  located  above  the  mine. 


Figure  2:  ALE  (left)  and  SPH  (right)  finite  element  models 


The  *Mat _plastic  kinematic  material  model  was  used  to  model  the  plate  and  the  frame,  and  the 
corresponding  material  properties  are  given  in  Table  1.  *Mat_high_explosive_burn  material 
model  and  the  *JWL  equation  of  state  (EOS)  were  used  for  the  explosive,  and  their  corresponding 
properties  are  provided  in  Table  2  and  Table  3  respectively.  *Mat_soil_and Joam failure 
material  model  was  chosen  for  the  soil  and  *Mat_null  material  model  with  a  linear  equation  of 
state  was  used  for  air.  Soil  and  air  material  parameters  are  given  in  Table  4,  Table  5  and  Table  6. 


Table  1:  Material  properties  used  in  the  experimental  setups  models 


Armour 
steel  1 

Armour 
steel  2 

Weld* 

A1  5083-H131 

Density,  (p)  kg/m3 

7850 

7850 

7850 

2770 

Elastic  modulus,  (E)  GPa 

197.5 

197.5 

197.5 

70.33 

PoissonS  ratio,  (V) 

0.33 

0.33 

0.33 

0.33 

Yield  stress,  (cv)  GPa 

1.0 

1.25 

1.04 

0.322 

Tangent  modulus,  (ET)  GPa 

2.175 

4.751 

1.391 

0.34 

Failure  strain,  (6f)  mm/mm 

N/R 

N/R 

0.225 

N/R 

*  [Bouamoul  and  Gaudreault,  2006] 
N/R:  Not  required 


Table  2:  *Mat_high_explosive_burn  material  properties  [Dobratz  and  Crawford,  1985] 


Explosive 

Density,  (p)  kg/m3 

1630 

Detonation  velocity,  (D)  m/s 

6930 

Chapman-Jouget  pressure,  (PCj)  GPa 

21 
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Table  3:  *EOS_JWL  equation  of  state  parameters  [Dobratz  and  Crawford,  1985] 


JWL 

A,  GPa 

371.2 

B,  GPa 

3.231 

R1 

4.15 

R2 

0.95 

Omeg,  (w) 

0.30 

Internal  energy  density,  (E0)  GPa  -  m3/m3 

7.0 

Initial  relative  volume,  (Vo) 

1.0 

Table  4:  *Mat_soil_and  J'oam ^failure  material  parameters  [Martec  Limited,  2007] 


Soil 

Density,  (p)  kg/m3 

1852 

Shear  modulus,  (G)  GPa 

40.68 

Bulk  modulus,  (K)  GPa 

50 

A0,  Al,  A2,  PC 

0.0,  0.0,  1.190,  0.0 

VCR,  REF 

0.0,  0.0 

cl, (2,  c3,  c4,  c5 

0,  -0.01861,  -0.05616,  -0.08934,  -0.1276 

c6,  c7,  c8,  c9,  clO 

-0.1547,  -0.2131,  -0.2594,  -0.2633,  -1.00 

PI,  P2,  P3,  P4,  P5 

0.0,  4.58xl06,  1.5xl07,  2.92xl07,  5.92xl07 

P6,  P7,  P8,  P9,  P10 

9.81xl07,  2.894x10s,  6.507xl08,  6.95xl08,  3.753xl010 

Table  5:  *Mat_null  material  parameters 


Air  set  1 

Air  set  2 

Density,  (p)  kg/m3 

1.293 

1.293 

Pressure  cutoff,  (Pc)  Pa 

-100 

0 

Table  6:  *EOS_linear_polynomial  equation  of  state  parameters  for  the  air 
[Bouamoul  and  Nguyen-Dang,  2008] 


Linear 

polynomial 

Polynomial  equation  coefficient, 

CO,  Cl,  C2,  C3  and  C6 

0.0 

Polynomial  equation  coefficient, 

0.4 

C4  and  C5 

Initial  internal  energy,  E0  (MPa) 

0.25 

Initial  relative  volume,  V0 

1.0 
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2.3  Results:  comparison  between  ALE  and  SPH  methods 


Figure  3  shows  a  comparison  of  the  predicted  deformation  histories  between  the  two  numerical 
models  and  the  experimental  test  at  the  same  location  in  the  plate.  The  piezo  pins  were  located 
15.24  cm  (6  in)  from  the  centre  of  the  plate  in  the  experimental  test.  The  two  numerical  models 
predict  the  same  vertical-displacement  profile  and  are  in  agreement  with  the  experimental  data. 


Figure  3:  Comparison  between  the  two  FE  models  and  the  experimental  results 

The  final  deformation  of  the  plate  was  measured  after  the  trial  and  a  maximum  deformation  of 
298  mm  was  measured.  This  residual  experimental  deformation  was  compared  with  the  output  of 
the  two  models  and  they  are  given  in  Table  7.  The  SPF1  is  the  approach  that  is  in  best  accordance 
with  the  final  experimental  deformation. 


Table  7:  Comparison  of  centre-plate  deformation  results 


Deformation  (mm) 

Relative  error 

Experimental 

298 

- 

SPH 

313 

5% 

ALE 

327 

9.7% 

2.4  Summary 

In  this  simple  scenario,  two  methods,  ALE  and  SPH,  were  used  to  compute  the  final  total 
deformation  of  a  plate.  The  same  profile  for  the  vertical  displacement  was  obtained  for  both 
approaches  and  was  in  agreement  with  experimental  data.  However,  SPH  is  the  approach  that 
was  in  best  accordance  with  the  final  experimental  deformation. 
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3  Second  scenario:  typical  LAV  mock-up 


3.1  Experimental  setup 

The  LAV  mock-up  is  a  testing  facility  developed  by  DRDC  Valcartier  to  study  the  reaction  of  the 
LAV  structure  and  occupants  when  subjected  to  side  and  under  belly  blast  attacks.  Figure  4 
shows  the  mock-up  that  was  developed  to  conduct  full-scale  trials  and  to  evaluate  the 
effectiveness  of  systems  such  as  structural  hull  reinforcements,  add-on  armour  systems,  seating 
and  stowage  systems.  The  mock-up  represents  a  2-m  long  section  of  the  rear  of  a  typical 
Canadian  LAV,  and  for  this  testing  configuration  the  total  mass  was  9600  kg.  It  consists  of  rigid 
frames  at  the  front  and  aft  of  the  mock-up,  on  which  is  bolted  a  sub-assembly  undergoing 
destructive  testing  and  representing  the  area  of  interest.  For  side  blast  trials,  the  sacrificial  target 
is  a  welded  assembly  of  the  sidewall,  sponson  top  panel,  sponson  wall  and  a  section  of  the  floor. 
The  mock-up  is  equipped  with  wheels  to  provide  representative  geometry  and  masking. 


Figure  4:  DRDC  Valcartier  LAV  mock-up  testing  facility 


The  mock-up  was  instrumented  with  five  3-axis  accelerometers  mounted  on  the  interior  walls  (at 
the  centre  of  the  sidewall,  at  the  centre  of  the  sponson  top,  on  the  comer  of  the  sponson,  at  the 
centre  of  the  sponson  wall  and  at  the  base  of  the  sponson  wall)  to  monitor  the  local  reaction  of  the 
structure.  The  accelerations  obtained  experimentally  were  integrated  to  generate  the  velocity 
history  curves,  which  were  compared  to  numerical  results.  Also,  a  laser  displacement  sensor  was 
aimed  at  the  comer  of  the  sponson.  Internal  high  speed  imagery  monitored  the  deformation  of  the 
wall  and  four  3 -axis  accelerometers  were  mounted  on  the  roof  to  measure  global  motion  of  the 
test  rig.  Wall  reflected  pressure  measurements  were  taken  inside  the  test  rig.  External  incident 
pressure  measurements  were  taken  at  10  m  and  20  m  from  the  charge.  Figure  5  shows  pictures  of 
the  internal  instrumentation  [ Toussaint  and  Durocher,  2008]. 
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Figure  5:  Internal  instrumentation 


3.2  Finite  element  models 

Figure  6  shows  the  ALE  and  SPH  FE  models  developed  to  support  the  experiments  performed  on 
the  mock-up  [Toussaint  and  Durocher,  2008].  In  both  models,  the  sub-assembly  parts  (mock-up 
frame)  were  generated  with  2,750  shell  elements  and  the  wheels  with  1,078  shell  elements.  The 
material  model  used  for  the  mock-up  frame  was  *Mat _plastic  kinematic.  Two  types  of  armour 
steel  were  used  for  these  parts  and  their  properties  are  given  in  Table  1.  The 
*Mat_high_expIosive_burn  material  model  and  the  *EOS_JWL  equation  of  state  were  used  for 
the  explosive,  and  their  corresponding  properties  are  given  in  Table  2  and  Table  3  respectively. 
The  wheels  were  modelled  for  masking  purposes.  The  material  model  *Mat_blatz-ko_rubber  was 
used  to  model  the  wheels,  and  Table  8  provides  their  corresponding  parameters.  In  addition,  a 
rigid  wall  condition  was  prescribed  at  the  ground  level  to  allow  reflection  of  the  ALE  air  blast  and 
the  SPH  particles  on  the  ground.  Finally,  the  upper  part  of  the  sidewall,  the  symmetry  plane  of 
the  floor  and  both  sides  of  the  mock-up  were  constrained  in  their  6  degrees  of  freedom. 


Figure  6:  ALE  (left)  and  SPH  (right)  finite  element  models 
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Table  8:  *Mat_Blatz-ko_rubber  properties  [Bouamoul  and  Nandlall,  2006] 


Rubber 

Density,  (p)  kg/m3 

1810 

Shear  modulus,  (G)  MPa 

2.53 

PoissonS  ratio,  v 

0.463 

Table  9:  *EOS_ideal_gas  equation  of  state  parameters  for  the  air  [Martec  Limited,  2007] 


Ideal  gas 

Nominant  constant-volume 
specific  heat  coefficient,  (Cv0) 

718 

Nominant  constant-pressure 
specific  heat  coefficient,  (Cp0) 

1005 

Initial  temperature,  (T0) 

270 

Initial  relative  volume,  (Vo) 

1 

For  the  ALE  model,  the  air  was  modelled  with  2,833,572  hexahedron  ALE  elements.  The 
dimensions  of  these  elements  were  0.02  m  length,  0.02  m  width  and  0.02  m  depth.  The  model 
was  large  enough  to  enclose  the  explosive  and  the  mock-up.  The  *Mat_nuII  and  the 
*EOS_ideal_gas  equation  of  state  were  used  to  model  the  air.  These  properties  are  summarized 
in  Table  5  and  Table  9  respectively.  The  explosive  geometry  was  defined  in  the  air  mesh  using 
the  *Initial_volume Jr action _geometry  feature  and  was  modelled  with  3,375  hexahedron  ALE 
elements.  A  ^Constrained _lagrange_in_solid  type  contact  was  defined  between  the  ALE 
elements  and  the  Lagrangian  structure.  Finally,  boundary  non-reflecting  conditions  were  used  in 
each  side  and  for  the  top  of  the  ALE  air  box  to  avoid  blast  wave  reflection  from  each  side  of  the 
ALE  box. 

For  the  SPFI  model,  the  explosive  was  modelled  with  215,940  SPF1  elements.  To  limit  the 
number  of  particles  to  be  traced  by  the  simulation  and  thus  enhance  the  calculation  time,  a  box 
envelope  was  specified  using  the  *C ontrol  sph  card.  In  this  way,  only  the  interactions  between 
particles  inside  this  box  were  considered.  Finally,  two  types  of  contact  were  used: 
*  Contact  automatic  nodes  between  SPH  particles  and  the  mock-up  structure  and 
*Contact_automatic_surface  between  all  the  Lagrangian  entities. 

3.2.1  ALE  method 

3. 2. 1.1  Mock-up  mesh  sensitivity  analysis 

A  mesh  sensitivity  analysis  of  the  mock-up  was  done  using  the  ALE  method.  In  the  first  mock- 
up  mesh,  the  sub-assemblies  (sidewall,  sponson  top  and  sponson  wall)  and  the  wheels  were 
modelled  with  2,750  (0.050  m  length  on  each  side  of  the  hexahedron  element)  and  1,078  shell 
Belytschko-Tsay  elements  respectively.  The  second  mock-up  mesh  contained  9,800  elements 
(0.025  m  length  on  each  side  of  the  hexahedron  element)  whereas  the  wheels  were  meshed  with 
1078  shell  Belytschko-Tsay  elements.  These  two  meshes  are  shown  in  Figure  7. 
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Figure  7:  Mesh  1  (left)  and  mesh  2  (right)  of  the  mock-up 


The  velocities  obtained  for  the  sidewall,  sponson  top  and  sponson  wall  with  mesh  1  and  mesh  2 
were  compared  to  the  experimental  measurements1.  Figure  8  shows  the  sponson  wall  horizontal 
velocity  (y-axis)  according  to  the  axis  orientation  shown  in  Figure  7.  Figure  9  shows  the  sponson 
top  vertical  velocity  (z-axis)  and  Figure  10  shows  the  sidewall  horizontal  velocity  (y-axis). 
Because  the  time  of  arrival  of  the  blast  was  not  of  interest  but  only  the  duration  and  the  maximum 
velocity,  the  beginning  of  the  velocity  curves  were  aligned  together  to  make  their  comparison 
easier.  This  alignment  was  also  done  for  all  the  following  figures  of  this  report. 


Figure  8:  Comparison  of  the  velocity  at  the  centre  of  the  sponson  wall  between  the  two  mock-up 

meshes  and  the  experimental  results 


1  It  is  important  to  mention  that  the  experimental  data  were  obtained  from  only  one  field  test  and  therefore, 
more  experimental  tests  may  be  needed  to  average  the  experimental  velocity  profile. 
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Figure  9:  Comparison  of  the  velocity  at  the  centre  of  the  sponson  top  between  the  two  mock-up 

meshes  and  the  experimental  results 


Figure  10:  Comparison  of  the  velocity  at  the  centre  of  the  sidewall  between  the  two  mock-up 

meshes  and  the  experimental  results 

As  seen  in  Figures  8,  9  and  10,  the  use  of  coarse  and  fine  meshes  gave  similar  velocity  profiles. 
However,  in  Figure  10,  the  numerical  peak  velocity  was  approximately  30  percent  less  than  the 
experimental  one.  A  combination  of  many  parameters  may  explain  this  difference:  the  use  of  the 
keyword  ^Boundary  non  reflecting  on  the  ALE  elements  boundary  instead  of  modelling  an 
infinite  domain,  the  air  material  model  and  equation  of  state  and  the  shadowing  effect  from  the 
wheels.  Also,  the  nodes  located  on  the  mock-up  boundaries  were  all  fixed  in  space,  which  did  not 
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allow  any  global  motion  of  the  mock-up.  This  was  not  the  case  in  the  experimental  trial,  where 
the  mock-up  was  free  to  move  with  the  blast. 

The  coarse  mock-up  mesh  was  used  in  the  subsequent  simulations  since  both  meshes  gave  similar 
results.  The  following  section  compares  two  techniques  used  to  simulate  the  ALE  air  blast. 

3. 2. 1.2  ALE  moving  mesh  box  study 

To  enhance  the  calculation  time  it  takes  to  run  the  mock-up  ALE  (Eulerian)  simulation,  an  ALE 
model  with  a  moving  domain  reference  was  generated  and  compared  to  a  fixed  domain 
(Eulerian).  The  difference  between  using  fixed  and  moving  domain  references  is  that  fewer  air 
elements  need  to  be  meshed  in  a  moving  domain  reference.  This  could  be  useful  when  the 
standoff  distance  between  the  explosive  and  the  structure  is  large,  which  means  that  more  air 
elements  are  required  to  enclose  all  parts  in  the  problem. 

The  FE  model  with  a  moving  domain  consisted  of  a  smaller  ALE  box  (896,700  hexahedron  ALE 
elements)  centred  initially  at  the  explosive  location.  This  box,  shown  in  Figure  1 1 ,  was  moving 
while  the  explosive  was  expanding  towards  the  structure.  The  * Ale  reference _sy stem _group  and 
*Ale_reference_system_curve  cards  were  used  to  allow  the  translation  of  this  box  towards  the 
mock-up.  The  velocity  of  the  box  was  adjusted  to  fit  the  expansion  of  the  explosive. 


Figure  11:  ALE  moving  mesh 

The  LS-DYNA  version  mpp971s  R2  7600.1224  and  four  32-bit  processors  were  used  to  perform 
the  ALE  moving  mesh  simulation,  while  the  ALE  fixed  mesh  (Eulerian  mesh)  simulation  was  run 
with  one  64b  processor  and  the  LS-DYNA  version  ls971s  R2  7600.1224.  We  expect  that  the 
numerical  results  obtained  using  these  two  different  platforms  (32b  versus  64b)  and  processors  (1 
versus  4)  may  be  slightly  different  but  still  not  significantly  so. 

To  verify  the  accuracy  of  the  results  generated  with  the  ALE  moving  domain  and  the  one 
generated  with  a  fixed  domain,  the  numerical  results  were  compared  to  the  experimental  ones. 
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Figure  12,  Figure  13  and  Figure  14  present  the  velocities  obtained  at  the  centre  of  the  sponson 
wall,  sponson  top  and  sidewall  of  the  mock-up  in  the  horizontal  (y),  vertical  (z)  and  horizontal  (y) 
direction  respectively  for  the  fixed  mesh,  the  moving  mesh  and  the  experimental  results. 


Figure  12:  Comparison  between  the  velocity  at  the  centre  of  the  sponson  wall  for  the  ALE  (fixed 

and  moving  mesh)  and  the  experimental  results 


Figure  13:  Comparison  between  the  velocity  at  the  centre  of  the  sponson  top  for  the  ALE  (fixed 
and  moving  mesh)  and  the  experimental  results 
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Figure  14:  Comparison  between  the  velocity  at  the  centre  of  the  sidewall  for  the  ALE  (fixed  and 

moving  mesh)  and  the  experimental  results 

These  figures  show  that  the  ALE  moving  mesh  slightly  underpredicted  the  velocities  of  the  fixed 
mesh  for  the  sponson  top  and  sidewall  but  greatly  overpredicted  the  sponson  wall  velocities. 
Consequently,  an  ALE  fixed  mesh  was  used  for  the  future  calculations.  Annex  A  gives  a 
summary  of  the  different  parameters  used  in  the  ^Constrained  lagrange  in  solid  card  for  these 
calculations. 

The  next  section  presents  a  comparison  of  three  ALE  fixed  air  meshes. 

3. 2. 1.3  Air  mesh  sensitivity  analysis 

This  section  presents  a  comparison  between  three  different  air  mesh  sizes  (mesh  3,  mesh  4  and 
mesh  5)1 2  for  the  ALE  fixed  mesh.  The  fifth  mesh  is  the  finer  one  and  was  the  one  used 
previously  to  perform  the  mock-up  mesh  sensitivity  analysis  and  the  ALE  moving  mesh  box 
study.  The  dimension  of  the  air  box  was  the  same  in  all  three  meshes  but  the  number  of  air 
elements  was  different.  The  third  mesh  contained  355,691  □  1  point  integration  ALE  elements 
with  approximate  dimensions  of  0.04  m  length,  0.04  m  width  and  0.04  m  depth.  The  fourth  mesh 
contained  831,870  □  1  point  integration  ALE  elements  with  approximate  dimensions  of  0.03  m 
length,  0.03  m  width  and  0.03  m  depth,  and  the  fifth  mesh  contained  2,833,572  □  1  point 
integration  ALE  elements  with  approximate  dimensions  of  0.02  m  length,  0.02  m  width  and  0.02 
m  depth.  A  cubic  geometry  was  used  to  model  the  explosive. 

Figure  15,  Figure  16  and  Figure  17  present,  respectively,  a  comparison  between  the  velocities 
obtained  at  the  centre  of  the  sponson  wall,  sponson  top  and  sidewall  in  the  horizontal  (y),  vertical 
(z)  and  horizontal  (y)  direction  for  three  air  meshes  and  the  experimental  results. 


1  The  air  meshes  were  named  mesh  3,  mesh  4  and  mesh  5  to  distinguish  from  mesh  1  and  mesh  2  used  for 
the  mock-up  meshes. 
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Figure  15:  Comparison  between  the  velocity  at  the  centre  of  the  sponsor  wall  for  three  air- 

meshes  and  the  experimental  results 


Figure  16:  Comparison  between  the  velocity  at  the  centre  of  the  sponsor  top  for  three  air  meshes 

and  the  experimental  results 
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Figure  17:  Comparison  between  the  velocity  at  the  centre  of  the  sidewall  for  three  air  meshes  and 

the  experimental  results 

As  seen  in  Figure  15  and  Figure  16,  mesh  3  and  mesh  4  provided  similar  results  with  maximum 
peak  velocities  close  to  the  experimental  data  while  mesh  5  showed  a  maximum  peak  velocity 
higher  than  the  experimental  one.  Figure  17  showed  similar  results  for  the  three  meshes  with 
peak  velocities  slightly  lower  than  the  experimental  one.  In  general,  fine  mesh  gave  better 
results,  however  in  this  case  mesh  3  and  mesh  4  provided  results  that  better  fit  with  the 
experimental  data.  This  could  be  explained  by  the  fact  that  in  mesh  5  the  number  of  coupling 
points  between  the  ALE  and  Lagrangian  meshes  was  higher  than  the  one  in  mesh  3  and  mesh  4. 
According  to  the  LS-DYNA  Keyword  UserS  Manual  [2007],  CToo  many  coupling  points  can 
result  in  instability,  and  not  enough  can  result  in  leakage  □  Therefore,  mesh  5  should  not  be  used 
in  conjunction  with  the  mesh  1  of  the  mock-up.  Other  factors  that  were  previously  mentioned  in 
Section  3. 2. 1.1  might  have  had  an  influence  on  the  numerical  results. 

The  simulation  took  1  hr  23  min.  for  mesh  3,  6  hrs  49  min.  for  mesh  4  whereas  mesh  5  took  29 
hrs  28  min.  Therefore,  to  optimize  the  time  required  to  run  a  simulation  and  considering  that 
mesh  3  and  mesh  4  provided  similar  results,  the  coarse  air  mesh  (mesh  3)  was  used  in  the 
subsequent  analysis. 

The  next  section  presents  the  effect  of  the  explosive  geometry. 

3. 2. 1.4  Effect  of  the  explosive  geometry 

A  numerical  study  was  performed  to  determine  if  the  initial  explosive  geometry  would  have  had 
an  influence  on  the  reaction  of  the  mock-up.  It  was  expected  that  the  geometry  of  the  explosive 
should  have  an  influence  at  close  range,  but  not  at  long  range. 

Three  different  explosive  geometries  were  modelled:  cube,  sphere  and  cylinder  (L/D  =1).  To 
preserve  the  total  weight  of  the  explosive,  the  volume  of  the  explosive  was  used  to  compute  the 


18 


DRDC  Valcartier  TR  2010-326 


dimensions  of  these  three  shapes.  The  FE  analyses  were  performed  using  the  coarse  air  mesh 
(mesh  3).  Figure  18,  Figure  19  and  Figure  20  present  the  velocities  obtained  at  the  centre  of  the 
sponson  wall,  sponson  top  and  sidewall  of  the  mock-up  in  the  horizontal  (y),  vertical  (z)  and 
horizontal  (y)  direction  respectively  for  three  different  explosive  geometries  and  the  experimental 
results. 


Figure  18:  Comparison  between  the  velocity  at  the  centre  of  the  sponson  wall  for  three  different 
explosive  geometries  and  the  experimental  results 


Figure  19:  Comparison  between  the  velocity  at  the  centre  of  the  sponson  top  for  three  different 
explosive  geometries  and  the  experimental  results 
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Figure  20:  Comparison  between  the  velocity  at  the  centre  of  the  sidewall  for  three  different 
explosive  geometries  and  the  experimental  results 

From  Figure  18,  Figure  19  and  Figure  20,  the  geometry  of  the  explosive  did  not  influence  the 
response  of  the  mock-up.  The  use  of  the  box  geometry  resulted  in  slightly  higher  velocities  than 
with  the  cylinder  or  the  sphere.  If  the  charge  had  been  closer  to  the  structure,  the  influence  of  the 
geometry  would  have  probably  been  more  important.  The  cube  shape  was  used  for  the 
subsequent  FE  analyses. 

3. 2. 1.5  Effect  of  the  air  EOS  parameters  and  the  air  pressure  cut-off  value 

The  purpose  of  this  section  is  to  compare  the  results  obtained  using  the  *Mat_null  material  model 
using  two  pressure  cut-off  values,  -100  Pa  and  0  (properties  are  given  in  Table  5)  and  to  compare 
the  results  obtained  using  two  equation  of  states  (EOS):  *EOS linear polynomial  and 
*EOS_ideal_gas  (properties  are  given  in  Table  6  and  Table  9  respectively). 

In  the  LS-DYNA  Keyword  UserS  Manual  [2007],  the  *E OS_l in ear _po lynom ia  1  equation  of  state 
is  defined  by  equation  (1). 

P  =  C0  +  C>  +  C0jU~  +  (73//3  +  (C4  +  C5ju  +  C6jU  (1) 

where  ja.=(p/po)- 1 ;  p  and  p0  are  the  density  and  the  reference  density  respectively;  C4  and  C5  are 
constants  to  be  defined  and  E  is  the  internal  energy  defined  in  unit  of  pressure.  If  pcO  then  Cfp 
and  C6/r  are  set  to  zero.  For  an  ideal  gas,  C0,  Ch  C2,  C3,  C6  are  set  to  zero  and  equation  1 
becomes: 


P  =  {C4+C5Ju)E 


(2) 
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By  introducing  the  value  of  po,  C4  and  C5  (given  in  Table  6),  equation  (2)  in  the  SI  system 
becomes: 

P  =  0.3094pE  (3) 

At  time  equals  zero,  the  initial  pressure  is  computed  by  introducing  the  reference  density  (given 
in  Table  5)  and  the  initial  internal  energy  (given  in  Table  6),  and  equals  100.014  kPa. 

In  the  LS-DYNA  Keyword  UserB  Manual  (2007),  the  *EOS_ideal_gas  equation  of  state  is 
defined  by  equation  (4). 

P  =  p(Cp-C,)t  (4) 

where  p  is  the  reference  density  (kg/m3);  Cp  and  Cv  are  the  specific  heat  capacity  at  constant 
pressure  and  volume  respectively  (J/(kg  K))  and  T  is  the  temperature  (K).  Assuming  that  Cp  and 
Cv  are  constants,  Cp  =  Cpo  and  Cv  =  Cv0  and  by  introducing  these  values  (given  in  Table  9)  in 
equation  (4),  equation  (4)  in  the  SI  system  becomes: 

P  =  287 pT  (5) 

At  time  equals  zero,  the  initial  pressure  is  computed  by  introducing  the  reference  density  (given 
in  Table  5)  and  the  initial  temperature  (given  in  Table  9),  and  equals  100.195  kPa. 

The  difference  between  the  initial  pressure  computed  from  equation  (3)  and  equation  (5)  is  equal 
to  0.181  kPa  which  is  negligible  (0.18%  relative  difference  between  the  two  pressures). 
Therefore,  similar  results  are  expected  when  choosing  one  equation  over  another  in  the  numerical 
model. 

Table  10  shows  a  summary  of  the  three  simulations  parameters.  For  the  first  case,  the 
*EOS  ideal _gas  equation  of  state  was  used  and  in  the  *Mat_null  material  model  the  pressure  cut¬ 
off  value  was  set  to  -100  Pa.  In  the  second  case,  the  same  EOS  was  used,  but  the  pressure  cut-off 
value  was  set  to  zero  and  for  the  third  case,  the  *EOS_linear_poIynomial  equation  of  state  was 
used  and  the  pressure  cut-off  value  was  set  to  zero. 


Table  10:  EOS  and  pressure  cut-off  value  for  each  case 


EOS 

Pressure  cut-off  value  (Pa) 

Case  1 

EOS_ideal_gas 

-100 

Case  2 

E  OS  idea  l_gas 

0 

Case  3 

EOS  linear  polynomial 

0 
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Figure  21,  Figure  22  and  Figure  23  present  the  velocities  obtained  at  the  centre  of  the  sponson 
wall,  sponson  top  and  sidewall  of  the  mock-up  in  the  horizontal  (y),  vertical  (z)  and  horizontal  (y) 
direction  respectively  for  the  three  cases  and  the  experimental  results. 


Figure  21:  Comparison  between  the  velocity  at  the  centre  of  the  sponson  wall  for  different  air  set 

parameters  and  the  experimental  results 


Figure  22:  Comparison  between  the  velocity  at  the  centre  of  the  sponson  top  for  different  air  set 

parameters  and  the  experimental  results 
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Figure  23:  Comparison  between  the  velocity  at  the  centre  of  the  sidewall  for  different  air  set 

parameters  and  the  experimental  results 

Figure  21,  Figure  22  and  Figure  23  show  that  the  pressure  cut-off  value  did  not  influence  the 
results,  i.e.  in  all  figures,  case  1  and  case  2  gave  exactly  the  same  velocity  profiles.  Flowever, 
contrary  to  expectations,  the  results  obtained  with  both  EOS  were  different.  The 
*EOS  linear _polynomial  equation  of  state  predicted  higher  velocities  than  the  *EOS_ideal_gas 
did.  According  to  the  LS-DYNA  Keyword  UserS  Manual  [2007]  and  for  the  Eulerian/ALE 
models,  the  *EOS  ideal _gas  equation  of  state  preserves  the  adiabatic  state  during  the  advection 
and  therefore,  the  internal  energy  is  perfectly  conserved.  This  might  explain  the  difference 
obtained  with  both  EOS.  Since  we  donl  have  much  information  about  the  feature,  it  is 
recommended  to  run  an  ALE  problem  with  both  EOS. 

Overall,  the  pressure  cut-off  value  did  not  influence  the  results.  Since  the  FE  model  with  the 
*EOS_linear_polynomial  equation  of  state  is  the  one  that  provided  better  results  compared  with 
the  experimental  data,  it  was  the  one  that  was  chosen  in  the  subsequent  analyses. 

In  summary,  the  set-up  that  was  selected  for  the  ALE  numerical  model  was:  the  mock-up  mesh  1 
(coarse  mesh),  a  fixed  ALE  domain,  the  ALE  mesh  3  (coarse  mesh),  box  geometry  for  the 
explosive  and  the  * E OSJ incur  jx) lyn omia l  equation  of  state  with  a  pressure  cut-off  set  to  zero. 

3.2.2  SPH  method 

To  optimise  the  numerical  model  using  SPF1  particles,  two  studies  were  performed.  The  first 
compares  three  different  SPH  density  meshes  and  the  second  compares  the  results  obtained  with 
different  particle  approximation  theories. 
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3.2.2. 1  Effect  of  the  number  of  SPH  particles 


A  mesh  sensitivity  analysis  was  performed  using  3,375-SPH  (mesh  6),  35,937-SPH  (mesh  7)  and 
215,940-SPH  (mesh  8)  particles  to  study  the  structural  deformations  of  the  mock-up.  Figure  24, 
Figure  25  and  Figure  26  present  the  velocities  obtained  at  the  centre  of  the  sponson  wall,  sponson 
top  and  sidewall  of  the  mock-up  in  the  horizontal  (y),  vertical  (z)  and  horizontal  (y)  direction 
respectively  for  the  three  SPF1  meshes  and  the  experimental  results. 
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Figure  24:  Comparison  between  the  velocity  at  the  centre  of  the  sponson  wall  for  three  SPH 

meshes  and  the  experimental  results 


Figure  25:  Comparison  between  the  velocity  at  the  centre  of  the  sponson  top  wall  for  three  SPH 

meshes  and  the  experimental  results 
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Figure  26:  Comparison  between  the  velocity  at  the  centre  of  the  sidewall  for  three  SPH  meshes 

and  the  experimental  results 

As  seen  in  Figure  24,  25  and  26,  the  coarse  mesh  led  to  more  oscillations  of  the  velocity  profiles 
than  the  fine  mesh  did.  These  oscillations  may  be  the  results  of  the  dispersion  of  the  particles 
near  the  structure  and  the  local  effect  of  the  particles  on  the  structure.  Also,  the  particles  did  not 
bypass  the  wheels  like  a  fluid  would  do,  so  it  was  difficult  to  reproduce  fluid-like  behaviour  of 
the  blast  on  the  structure.  This  effect  was  also  observed  in  the  work  done  by  Toussaint  and 
Durocher  [2008].  Finally,  in  this  SPFI  model,  the  air  was  not  modelled. 

The  finer  mesh  (215,940  SPFI  particles)  is  the  one  that  provided  the  best  fit  with  the  experiment 
and  thus,  was  used  in  the  subsequent  analysis. 

3.2. 2. 2  Effect  of  the  particle  approximation  theory 

Seven  particle  approximation  theories  are  available  in  LS-DYNA  code.  They  are  listed  in 
Table  12.  All  seven  formulations  were  simulated.  Formulations  1,  3,  4  and  6  computed  infinite 
velocities,  therefore  the  velocity  profiles  were  only  given  for  formulations  0,  2  and  5.  Figure  27 
Figure  28  and  Figure  29  present  the  velocities  obtained  at  the  centre  of  the  sponson  wall,  sponson 
top  and  sidewall  of  the  mock-up  in  the  horizontal  (y),  vertical  (z)  and  horizontal  (y)  direction 
respectively  for  the  three  SPFI  formulations  and  the  experimental  results. 


DRDC  Valcartier  TR  2010-326 


25 


Table  11:  Particle  approximation  theory  [LS-DYNA  Keyword  Userid  Manual,  2007] 


Formulation 

number 

Particle  approximation  theory 

Simulation 

0  (default) 

Not  available 

yes 

1 

Renormalization  approximation 

no 

2 

Symmetric  formulation 

yes 

3 

Symmetric  renormalized 

4 

approximation 

Tensor  formulation 

no 

5 

Fluid  particle  approximation 

yes 

6 

Fluid  particle  with 
renormalization  approximation 

no 

Figure  27:  Comparison  between  the  velocity  at  the  centre  of  the  sponson  wall  for  formulation  2,  5 

and  the  default  and  the  experimental  results 


26 


DRDC  Valcartier  TR  2010-326 


60.00  -i 


40.00 


>. 

o 

O 

o 

> 


20.00 


0.00 


-20.00 


-40.00 


Time  (s) 


Figure  28:  Comparison  between  the  velocity  at  the  centre  of  the  sponson  top  for  formulation  2,  5 

and  the  default  and  the  experimental  results 


Figure  29:  Comparison  between  the  velocity  at  the  centre  of  the  sidewall  for  formulation  2,  5  and 

the  default  and  the  experimental  results 

As  shown  in  Figure  27,  formulations  0  and  2  gave  similar  peak  velocities  and  formulation  5 
provided  the  highest  peak  velocity,  but  all  three  formulations  computed  velocities  that  were  far 
from  the  experiment.  Some  factors  that  could  explain  these  results  were  listed  in  section  3.2.2. 1. 
In  Figure  28,  all  three  formulations  computed  similar  peak  velocities,  but  again,  they  were  far 
from  the  experimental  data.  In  Figure  29,  formulation  2  is  the  one  that  was  close  to  the 
experimental  data  and  thus  is  the  one  that  provided  the  best  fit  with  the  experimental  data. 

In  summary,  the  set-up  that  was  selected  for  the  SPH  numerical  model  was:  the  finer  mesh 
(215,940  SPH  particles)  and  formulation  2  for  the  particle  approximation  theory. 
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3.3  Results:  comparison  between  ALE  and  SPH  methods 


To  determine  which  method  ALE  or  SPH  was  more  suitable  to  model  the  effect  of  air  blast  on  a 
structure,  a  comparison  between  the  best  outcome  of  the  ALE  and  the  SPH  formulations  was 
done  in  this  section.  The  ALE  FE  analysis  was  performed  using  LS-DYNA  version  ls971s  R2 
7600.1224  and  one  processor,  whereas  the  SPH  analysis  was  performed  using  LS-DYNA  version 
mpp  971s  R2  7600.1224  and  four  processors.  The  ALE  analysis  was  run  on  only  one  processor 
because  of  the  size  of  the  model  (a  large  amount  of  memory  was  necessary  to  run  the  analysis). 
This  required  running  the  ALE  simulation  on  a  64b  processor  and  only  one  64b  processor  was 
available. 

For  both  analyses,  the  same  FE  mock-up  model  (coarse  mesh)  was  used  and  the  explosive  was 
modelled  as  a  box.  For  the  ALE  analysis,  the  ALE  coarse  mesh  was  used  (355,691  □  1  point 
integration  ALE  elements)  with  a  fixed  domain  (Eulerian  mesh).  A  pressure  cut-off  was  set  to 
zero  in  the  *Mat_null  air  material  model  and  the  *EOS_linear_polynomial  equation  of  state  was 
used.  For  the  SPH  analysis,  the  numerical  simulation  was  run  with  an  explosive  containing 
215,940  SPH  particles  in  conjunction  with  formulation  2  for  the  particle  approximation  theory. 

It  took  1  hr  23  min.  to  run  the  ALE  problem  (coarse  mesh)  compared  to  47  min.  for  the  SPH 
model  (fine  mesh). 


Figure  30,  Figure  31  and  Figure  32  show,  respectively,  the  velocity  in  the  horizontal  (y)  direction 
at  the  centre  of  the  sponson  wall,  the  velocity  in  the  vertical  (z)  direction  at  the  centre  of  the 
sponson  top  and  the  velocity  in  the  horizontal  (y)  direction  at  the  centre  of  the  sidewall  of  the 
ALE  and  SPH  numerical  model  results  compared  to  the  experimental  data. 


Figure  30:  Comparison  of  the  velocity  at  the  centre  of  the  sponson  wall  for  ALE,  SPH  and  the 

experimental  results 
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Figure  31:  Comparison  of  the  velocity  at  the  centre  of  the  sponson  top  for  ALE,  SPH  and  the 

experimental  results 


Figure  32:  Comparison  of  the  velocity  at  the  centre  of  the  sidewall  for  ALE,  SPH  and  the 

experimental  results 

The  SPH  model  overpredicted  the  sponson  wall  peak  velocity  and  underpredicted  the  sponson  top 
peak  velocity  compared  to  the  experimental  results.  The  sidewall  peak  velocity  was  similar  to  the 
experimental  results.  The  ALE  model  gave  similar  peak  velocities  for  the  sponson  wall  and 
sponson  top  but  underpredicted  the  sidewall  peak  velocity.  Overall,  the  velocities  obtained  with 
the  ALE  model  were  closer  to  the  experimental  velocities  than  those  obtained  with  the  SPH 
model.  Some  factors  could  explain  these  results.  In  the  SPH  meshless  formulation  when 
particles  are  far  from  each  other,  some  numerical  instability  may  occur.  This  effect  was 
emphasized  in  the  second  scenario  where  the  blast  was  far  from  the  structure.  Also,  shadowing  is 
not  accounted  for  when  using  SPH  formulation.  Also,  the  air  was  not  modelled  in  the  SPH 
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model.  For  all  these  reasons,  SPF1  is  probably  not  the  accurate  approach  to  use  when  the  air  blast 
is  far  from  the  structure. 


3.4  Summary 

In  the  mock-up  scenario,  different  numerical  analyses  were  performed  for  both  methods.  For  the 
ALE  model,  initially  a  mesh  sensitivity  analysis  of  the  mock-up  was  done  and  since  both  meshes 
were  in  agreement  with  the  experimental  data,  the  coarse  mesh  of  the  mock-up  was  the  one  that 
was  used  in  the  subsequent  analyses.  A  comparison  between  moving  and  fixed  ALE  air  meshes 
was  done  and  it  was  shown  that  for  the  current  application,  it  was  more  appropriate  to  use  an  ALE 
fixed  mesh  (Eulerian  mesh).  Also,  the  effect  of  air  element  size  was  studied  and  the  coarse  air 
mesh  (mesh  3)  was  chosen.  Then  the  effect  of  the  explosive  geometry  was  simulated,  and  it  was 
shown  that  for  this  particular  experimental  trial  and  with  the  actual  standoff  distance,  the 
geometry  of  the  explosive  did  not  significantly  influence  the  response  of  the  mock-up.  Finally, 
the  effect  of  air  material  properties,  especially  the  pressure  cut-off,  and  the  equation  of  state  were 
investigated.  It  was  demonstrated  that  for  the  ALE  model,  the  value  of  the  pressure  cut-off  in  the 
material  model  did  not  influence  the  results.  However,  different  results  were  obtained  with  the 
*EOS  linear polynomial  and  the  *EOS_ideal_gas  equations  of  state.  The 
*EOS  linear polynomial  was  chosen  since  it  provided  results  that  were  closer  to  the 
experimental  data.  An  SPH  mesh  sensitivity  analysis  was  performed  and  showed  that  a  fine  mesh 
(215,940-SPH  particles)  was  necessary  to  simulate  a  fluid-like  behaviour  and  to  avoid  local 
effects  (i.e.  when  particles  are  too  large,  they  impact  shell  elements  and  create  local  holes  instead 
of  creating  a  global  deformation  of  the  structure).  Also,  three  particle  approximation  theories 
were  investigated  and  formulation  2  is  the  one  that  provided  best  results. 

Factors  other  than  the  one  presented  in  this  study  could  have  influenced  the  numerical  results. 
For  example,  blast  shadowing  was  dependent  on  how  the  wheels  were  modelled  (dimensions  and 
properties)  and  thus  could  have  an  influence  on  the  numerical  results.  Also,  the  nodes  located  at 
the  outside  boundaries  of  the  mock-up  were  all  fixed  in  space,  which  did  not  allow  any  global 
motion  of  the  structure  (only  the  local  reaction  was  studied  in  this  work).  This  was  not  the  case  in 
the  experimental  trial,  where  the  mock-up  was  free  to  move  with  the  blast. 

Finally,  the  best  results  obtained  for  each  method  were  compared  and  it  was  shown  that  the  ALE 
formulation  reflected  the  expansion  of  the  air  blast  better  than  the  SPH  did.  Also,  since  only  one 
set  of  experimental  data  was  generated  so  far,  additional  tests  under  the  same  conditions  would  be 
required  to  assess  the  experimental  variation. 
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4  Conclusion 


One  major  issue  for  numerical  analysis  is  to  accurately  predict  the  mine  blast  effects  on  light 
armoured  vehicles  (LAV).  Therefore,  two  approaches  were  evaluated  in  this  report:  smoothed 
particle  hydrodynamics  (SPH)  and  arbitrary  Lagrangian  Eulerian  (ALE)  formulations  to  simulate 
a  mine  blast  on  a  typical  structure.  Two  scenarios  were  studied  using  these  two  methods. 

In  the  target  plate  scenario,  the  same  vertical-displacement  profile  was  obtained  for  both  ALE  and 
SPH  approaches  and  was  in  agreement  with  experimental  data.  However,  SPH  is  the  approach 
that  was  in  best  accordance  with  the  final  experimental  deformation. 

In  the  mock-up  scenario,  many  studies  were  performed  but  in  general  the  ALE  model  reflected 
the  evolution  of  blast  better  than  the  SPH  model  did. 

The  more  relevant  aspects  to  consider  when  using  the  ALE  method  to  model  an  air  blast  are  to 
use  a  fixed  domain  with  a  good  mesh  and  to  carefully  choose  the  air  material  properties  and  its 
corresponding  equation  of  state. 

In  general,  many  factors  may  influence  ALE  or  SPH  finite  element  results.  They  include:  the 
complexity  and  dimensions  of  the  structures  exposed  to  the  blast,  the  mesh,  the  material 
constitutive  models  and  their  corresponding  equation  of  state  and  the  boundary  conditions.  All 
these  parameters  need  to  be  investigated  in  more  detail  when  simulating  a  mine  blast.  Also,  since 
only  one  set  of  experimental  data  was  generated  so  far  for  the  mock-up,  additional  experimental 
tests  under  the  same  conditions  would  be  required  to  assess  the  experimental  variation. 

More  work  needs  to  be  done  in  order  to  determine  the  accurate  loading  to  be  used  for  landmine 
and  air  blast  scenarios.  Further  numerical  simulations  and  experiments  should  be  performed  to 
get  a  more  thorough  study  on  the  ALE  and  the  SPH  approaches.  In  particular,  it  would  be 
important  to  compare  not  only  the  local  reaction  of  the  finite  element  structure  with  experimental 
data  but  also  the  global  motion.  Also,  an  investigation  on  the  new  feature  available  in  LS-DYNA 
that  combines  the  *Load_blast_enhanced  card  and  the  ALE  method  could  be  relevant.  Finally,  it 
would  be  of  interest  to  test  the  ability  of  the  CHINOOK  code  to  model  landmine  scenarios. 
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Annex  A  ALE  parametric  study 


The  objective  was  to  find  the  best  combination  of  parameters  that  avoid  or  limit  the  leakage3  in 
the  mock-up  ALE  model.  The  different  parameters  studied  were:  NQUAD,  DIRECT,  PFAC, 
FRCMIN,  ILEAK  and  PLEAK.  The  definitions  of  these  parameters  are  given  in  the  LS-DYNA 
Keyword  UserH  Manual  [2007]  and  are  listed  below: 

NQUAD:  Number  of  coupling  points  distributed  over  each  coupled  Lagrangian  surface  segment. 
CTYPE:  Fluid-Structure  coupling  method. 

DIRECT :  Coupling  direction. 

PFAC:  Penalty  factor. 

FRCMIN:  Minimum  volume  fraction  of  a  coupled  ALE  multi-material  group  or  fluid  in  a  multi¬ 
material  ALE  element  to  activate  coupling. 

ILEAK:  Coupling  leakage  control  flag. 

IPLEAK:  Leakage  control  penalty  factor. 

Table  A-l  shows  the  different  parameter  combinations  and  the  yellow  row  gives  the  best  set  of 
parameters  that  was  used  in  the  ALE  numerical  results.  It  was  suggested  in  the  LS-DYNA 
Keyword  User H  Manual  [2007]  to  set  FRCMIN  between  0.1  and  0.3  in  hypervelocity  cases,  so 
FRCMIN  was  set  to  0. 1  since  it  turns  on  coupling  earlier  which  helps  preventing  leakage. 

Table  A-l:  ^Constrained _lagrange_in solid  parameters  study 


NQUAD 

CTYPE 

DIRECT 

PFAC 

FRCMIN 

ILEAK 

IPLEAK 

Leakage 

3 

4 

default 

0.2 

0.3 

default 

default 

a  lot 

3 

4 

3 

0.2 

0.3 

default 

default 

very  little 

3 

4 

3 

0.2 

default 

default 

default 

very  little 

3 

4 

default 

0.2 

0.1 

default 

default 

a  lot 

3 

4 

2 

0.2 

default 

default 

default 

some 

3 

4 

2 

0.2 

0.3 

default 

default 

some 

3 

4 

2 

0.2 

0.1 

default 

default 

some 

4 

4 

2 

0.2 

0.1 

default 

default 

a  lot 

4 

4 

2 

0.2 

0.3 

default 

default 

a  lot 

4 

4 

2 

0.2 

default 

default 

default 

a  lot 

4 

4 

2 

0.2 

0.1 

2 

0.2 

a  lot 

3 

4 

3 

0.2 

0.1 

2 

0.2 

N/A 

3  Leakage  appears  when  the  material  flow  passes  through  some  of  the  elements  of  the  Lagrangian  mesh. 
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List  of  symbols/acronyms 


ALE 

DND 

DRDC 

DRDKIM 

EEI 

EF 

FE 

hr  /  hrs 

IED 

LAV 

min 

RDDC 

R&D 

SPH 

VBL 

L/D 

P 

v 

G 

E 

Gy 

ETAN 

FS 

Vy 

vz 


Arbitrary  Lagrangian  Eulerian  /  Lagrangienne  Lul  ricnnc  arbitraire 

Department  of  National  Defence 

Defence  Research  and  Development  Canada 

Director  Research  and  Development  Knowledge  and  Information 
Management 

Engin  explosif  improvis  □ 

m  mcnts  finis 

Finite  element 

Hour  /  Hours 

Improvised  explosive  device 
Light  armoured  vehicle 
Minutes 

Recherche  et  d  Vcloppcmcnt  pour  la  defense  Canada 

Research  and  Development  /  Recherche  et  d  vcloppcmcnt 

Smoothed  particle  hydrodynamics  /  Hydrodynamique  des  particules  lisses 

VEhicule  blind  Dither 

Length  /  Diameter 

Mass  density 

Poisson  B  ratio 

Shear  modulus 

Young  S  modulus 

Yield  stress 

Tangent  modulus 

Failure  strain 

Velocity  in  the  horizontal  direction 
Velocity  in  the  vertical  direction 
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